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We propose a quasi-particle model for the thermodynamic description of the gluon plasma which takes into account non-abelian 
characteristics of the gluonic field. This is accomplished utilizing massive non-linear plane wave solutions of the classical equations 
of motion with a variable mass parameter, reflecting the scale invariance of the Yang-Mills Lagrangian. For the statistical description 
of the gluon plasma we interpret these non-linear waves as quasi-particles with a temperature dependent mass distribution. Quasi- 
Gaussian distributions with a common variance but different temperature dependent mean masses for the longitudinal and transverse 
modes are employed. We use recent Lattice results to fix the mean transverse and longitudinal masses while the variance is fitted 
to the equation of state of pure S U(3) on the Lattice. Thus, our model succeeds to obtain both a consistent description of the gluon 
plasma energy density as well as a correct behaviour of the mass parameters near the critical point. 
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1. Introduction 

Strongly interacting non-abelian gauge theories are described 
by S U(N) algebras. These theories possess a strong coupling in 
the low temperature regime prohibiting perturbative treatment. 
In addition, the associated degrees of freedom occur exclusively 
in a confined phase. In higher temperatures, near and above the 
critical point of S U{N) gauge theories, it is expected that the in¬ 
teraction of fermions and gauge fields, namely quarks and glu¬ 
ons in the case of S U(3) color (QCD), is significantly weaker, 
leading to the deconfined phase known as quark gluon plasma 
(QGP)!!!]. The thermodynamical treatment of QGP has led to 
the introduction of quasi-particle models (QPMs) Hi, pri¬ 
marily aiming to explain the QGP equation of state as obtained 
from Lattice gauge theory simulations of QCD at finite tem¬ 
perature. In these models and at temperatures higher than the 
critical value it was assumed that the thermodynamics of a sys¬ 
tem of interacting massless gluons may be approached by an 
ideal gas of massive noninteracting gluons. However, the above 
models |01 failed to explain the most recent Lattice results j^]. 


Due to asymptotic freedom, the interaction of quarks and 
gluons is expected to be very weak at extremely high tempera¬ 
tures. Thus, in these thermodynamic conditions the picture of 
an ideal gas of (almost) non-interacting particles for the quark- 
gluon system should provide a good approximation. Neverthe¬ 
less, as supported by the results of the experiments at the Rela¬ 
tivistic Heavy Ion Collider (RHIC) |@1, in the neighbourhood of 
the critical temperature Tc, associated with the transition from 
hadronic matter to QGP, the interaction is strong and the quark- 
gluon system is far from the ideal gas scenario 00] sharing 
features of a perfect fluid |@1. This holds in particular also 
for the gluon field alone where accurate Lattice results fin 
demonstrate that the gluon system remains far from the ideal 
behaviour even for temperatures 5 times larger than Tc- 


A useful and common strategy is to restrict the analysis to 
the gluonic sector considering the emergence of a (non-ideal) 
gluon plasma above the associated critical temperature. To cap- 
ture this non-ideal behaviour the QPMs ||iS[Iil[l2, isi 14, 15 


II [BUS ms introduce temperature dependent parame- 
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ters which are suitably adjusted in order to fit the existing Lat¬ 
tice results. A basic assumption of these models is the pres¬ 
ence of a temperature dependent mass for the gluons, a property 
which may lead to thermodynamic inconsistencies 12 111 within 
the Landau statistical approach 12^ . This is due to the fact 
that the temperature dependent mass becomes a thermodynamic 
quantity affecting the usual relations connecting pressure with 
energy density. This inconsistency may be healed by introduc¬ 
ing an appropriate constraint in models involving vacuum en¬ 
ergy B{T) 12111 . However, a much more natural way to over¬ 
come this issue is to use the Pathria approach starting the 
calculations of thermodynamic quantities and equation of state 
from the energy density instead of the pressure 1241] . The pres¬ 
sure is obtained via the integration of a fundamental thermody¬ 
namic relation, taking into account the temperature dependence 


of the gluon mass 11241] . The latter is determined by fitting the 
Lattice results for the equation of state of pure S U{3) at finite 
temperature ||3l . With the suitable temperature dependence for 
the gluon mass, the description of Lattice data with such an im¬ 
proved QPM turns out to be quite satisfactory. However, despite 
of being free from thermodynamic inconsistencies QPMs still 
include controversial issues from the physical point of view. 
The use of massive gluons and their treatment as free particles 
is not fully justified. To resolve this issue, a recent work lE^ 
assumed that the gluon mass emerges through the propagation 
of gluons in a plasma environment as a collective effect and is 
furthermore related to the associated plasma frequency. Lat¬ 
tice results for the equation of state are sufficiently reproduced 
with the use of a single temperature dependent parameter. Re¬ 
garding the origin of gluon mass, a consistent interpretation is 
possible in terms of classical solutions of the corrensponding 
gauge theory equations of motion (e.g see ll^ and references 
therein). In particular, 02711 derives a class of nonlinear plane 
wave S U(2) solutions obeying a massive relativistic dispersion 
relation. The mass parameter is free to vary as a consequence of 
the scale invariance of the Yang-Mills Lagrangian. Despite the 
fact that such a scenario is based on fundamental properties of 
Yang-Mills dynamics it has not been associated with the gluon 
mass in QPMs up to now. 

A novel aspect in the framework of QPMs emerged by re¬ 
cent Lattice S U{3) calculations of the temperature dependence 
of the dynamically generated gluon mass by estimating the in¬ 
verse gluon propagator in the infra-red limit ll2^ . These cal¬ 
culations obtained the gluon mass for a temperature regime just 
above the critical point, providing further constraints in the phe¬ 
nomenology of QPMs. There are some important consequences 
of the Lattice results: firstly the dynamical masses of the trans¬ 
verse and longitudinal gluon degrees of freedom differ at low 
temperatures. Secondly both masses (longitudinal and trans¬ 
verse) behave smoothly as a function of the temperature just 
beyond the critical point. Thus, the stiff increase of quasipar¬ 
ticle masses as the critical region is approached from above, is 
incompatible with the Lattice lE^ . A solution to this contro¬ 
versy has been proposed in the framework of QPMs through a 
Polyakov loop coupling to the quasiparicles 13011 . 

In the present work we develop a QPM which takes into ac¬ 
count the most recent Lattice results for the temperature de¬ 


pendence of the dynamical gluon mass 112911 and at the same 
time is in consistency with the older Lattice calculations of the 
equation of state ||3l. The proposed model suggests that the 
non-linear plane wave solutions of the equations of motion cor¬ 
respond to quasi-particles with variable mass. We develop such 
a scenario based on a subset of classical solutions for the gauge 
field, namely those which originate from the S t/(2)-sector con¬ 
tained in SU{3), assuming that the main characteristics of its 
non-abelian character are captured by this class. We demon¬ 
strate that the non-abelian character of the gluons introduces 
significant changes in their thermodynamical treatment which 
are taken into account in the proposed QPM. In contrast to 
other statistical models which use glueballs with a discrete mass 
spectrum Bill here we assume a continuously varying gluon 
mass characterized by a specific probability density. To take 
into account the difference between transverse and longitudi¬ 
nal masses, as calculated from the Lattice, we use two different 
mass distributions for transverse and longitudinal gluonic de¬ 
grees of freedom respectively. Furthermore, since the gluon 
mass distributions are unknown we choose truncated Gaussians 
(negative mass values are excluded) in each case. The Lattice 
results (extrapolating to high temperatures when necessary) ||3 
are employed to fix the temperature depended mean transverse 
and longitudinal masses. To reduce the number of free param¬ 
eters we further assume that the variation is the same in both 
quasi-Gaussian distributions. Thus, the temperature depended 
variance is the only free parameter in our model which is deter¬ 
mined via a fit to the Lattice results for the S U(3) equation of 
state ||3l. The main success of the proposed non-abelian quasi¬ 
particle model (NAQPM) is the very good description of two 
different Lattice results for the gluon plasma using a single free 
parameter. 

The paper is organized as follows. In Section II we present 
the non-linear plane wave solutions of the SU(2) Yang-Mills 
theory and we reveal their properties which are relevant for the 
subsequent formulation of the NAQPM. In section III we in¬ 
troduce the NAQPM for the gluon plasma and the correspond¬ 
ing statistical treatment. In section IV we discuss our results 
concerning the dependence of the parameters of the proposed 
model on temperature as well as their compatibility with exist¬ 
ing Lattice results. Finally, in section V we give our concluding 
remarks. 


2. Non linear plane waves in Sf7(2) Yang-Mills theory 

A simplified description of the gluon field, capturing the ba¬ 
sic phenomenological characteristics, can be obtained by the 
classical S U{2) Yang-Mills theory. Since S U{2) is a subgroup 
of S U{3), we will restrict our analysis to this case which is eas¬ 
ier to handle and at the end of this section we will give some ar¬ 
guments supporting that our treatment can be transferred to the 
more realistic description with colored S U(3) Yang-Mills. Ne¬ 
glecting fermionic (matter) degrees of freedom the Lagrangian 
density of this model is written as: 

£ = (d 


2 




where is the antisymmetric field tensor - 

dvA“ - g ■ eabcA^Ay of the gauge field The corresponding 
classical equations of motion for the gauge field are; 

d^F^y+g-eabc-A^^F^.^O ( 2 ) 

The equations ^ permit plane wave solutions of the form: 

Af; = i-e^OCwf-^-x) (3) 

g 

where <1) is a scalar function of the plane wave phase 

^ - oj ■ t - k ■ X (4) 


Solutions of the form of Eq. (|7]i are also found in the ((>^ scalar 
field theory forming a complete (non-orthogonal) basis 11321] for 
this system. Since there is a mapping of Yang-Mills S U(N) to 
the scalar (f)^ theory 11331] one may argue that the solutions in 
Eq. (I7]l provide also a non-orthogonal basis for SU{3) Yang- 
Mills gauge theory. The question is if and how the non-linear 
plane waves (l7]i can capture the main features of gluon plasma 
thermodynamics. Undoubtedly they will influence the counting 
of gluon microstates. The period of the non-linear plane waves 
is not 2n but V — 5.244, thus the number of such stationary 
states fitting in a fixed volume V is greater that the correspond¬ 
ing number for linear plane waves by a factor . 


and e" (for a = 1,2,3 and yu = 0,1,2,3) is given by 
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The momentum four-vector = (w, k) satisfies the dispersion 
relation —P' + iP for an arbitrary mass parameter m. Note 
that the three columns {a - 1,2,3) of the above matrix © are 
non-other than the three orthonormal spacelike vectors or¬ 
thogonal also to the timelike four-momentum vector k^j. As 
such, the solution obeys automatically the Lorentz gauge con¬ 
dition ^yjA^" = 0 since for plane waves it becomes equivalent to 
the transversality condition k^jA^" = 0. 

The solution © is most conveniently derived by the authors 
of ll^ on the proper time frame k^ = (m, 0) with the gauge 
fixing condition A" = 0. On that frame, Aj = Aj = = 

(m/g)<l>(^) with all other components equal to zero and d)(^) 
satisfying the equation; 


H- 2 ■ = 0 . (6) 


The above equation possesses solutions of the form: 

<I)(^) = 5n[^;-l]=cn[V2^;i] (7) 

where sn[^; k^](cn[^; k^]) is the Jacobi elliptic sine (cosine) 
function with elliptic modulus k lE^ . Thus, the solutions 
© describe periodic (anharmonic) plane waves with period 
V - ^K{k)l yfl where K{k) is the complete elliptic integral of 
the first kind. Eor = 1 /2, the period numerically becomes 
F - 5.244. A Lorentz boost with parameters j - calm, jS - 
kjcj leads immediately to the general form ©. The origin of 
the mass of the gluon field m is traced in the non-linear terms 
in the Lagrangian © due to the gluon self-interaction. Eurther- 
more, m is free to take any positive value, reflecting the scale 
invariance of Eq. ©. Thus, although scale-free, the Yang-Mills 
classical solutions depend on an arbitrary mass-scale m due to 
the nonlinearity of the theory. 


* Note that holds. 


Let us now proceed with the formulation of the NAQPM as 
discussed in the introduction. It is useful to list the main as¬ 
sumptions of the model; 

• The microstates of the gluon field at thermal equilibrium 
consist from non-linear plane wave solutions of the S U(2) 
classical equations of motion with period F = 5.244. 
Such classical field configurations in quark-gluon plasma 
studies have been used before ifjdl] . Within our treatment 
the classical S U{2) non-linear plane waves correspond to 
massive quasi-particles with a variable mass m. Since the 
S U(3) algebra contains three S U{2) subspaces, there ex¬ 
ist three different ways to embed the above S U{2) solution 
in S U{3). The vanisihing trace condition reduces by one 
the diagonal degrees of freedom, thus we effectively iden¬ 
tify eight degrees of freedom, in consistency with the full 
S U(3) case. In a more precise sense, the NAQPM does not 
consider the entire S U(3) solution space but rather the lo¬ 
cally isomorphic case of S U{2)xS U(2)xS I/(2)/t/(l) al¬ 
gebra. Nevertheless, this space contains the basic ingredi¬ 
ent of the non-abelian character, namely the non-linearity 
due to the gauge field self-interaction. A recent publica¬ 
tion has found a larger class of S U(3) plane wave so¬ 
lutions which could in principle be included in the present 
model. However, these solutions possess an infinite count¬ 
able set of periods in contrast to the solutions in the S U{2) 
subspace which have fixed period F - 5.244. Thus, their 
inclusion is highly non-trivial and goes beyond the scope 
of the present work. 

• Lattice calculations provide evidence that the transverse 
and longitudinal gluonic degrees of freedom acquire tem¬ 
perature depended masses which differ at low tempera¬ 
tures and approach each other as the temperature increases 
]|2^ . Eor sufficiently high temperatures these masses at¬ 
tain a common asymptotic linear dependence on the tem¬ 
perature in accordance with perturbation theory predic¬ 
tions i^ . In the NAQPM this information is incorpo¬ 
rated assuming quasi-Gaussian (with restriction to posi¬ 
tive values only) distributions for the gluon mass m which 
have well determined but different, temperature depen¬ 
dent, mean values gL,r{T) ~ {m)tr and nioiT) ~ {m)io 
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for the transverse (tr) and longitudinal (lo) gluonic de¬ 
grees of freedom respectively. In each case the average 
is meant over the microstates of the corresponding gluon 
components (transverse or longitudinal) with temperature 
depended weights. In such a description the associated 

variances cr,r ~ yj{m^),r-l4r and crio ~ yj{m^)io - 
depend also on temperature and present mass related re¬ 
sponse functions of the gluonic system. 


Within the NAQPM the temperature depended parameters 
jJir and yU/o are hxed by the recent lattice results ll29n on 
the gluon propagator. Since we need to calculate thermo¬ 
dynamic properties of the gluon system also for tempera¬ 
tures beyond those investigated in lE^ . we will appropri¬ 
ately extrapolate and /r/o. This extrapolation takes into 
account that in the high temperature regime both quanti¬ 
ties approach each other and attain asymptotically a linear 
temperature dependence as dictated by perturbation theory 
1^ . The variances CTtr and crio are the free parameters of 
the model which are determined in order to fit the Lattice 
results on the equation of state of the S U(3) Yang-Mills 
theory Etl. To reduce the number of free parameters we 
further assume that cr,r{T) - crio(T) = cr(T). We con- 
hrmed numerically that allowing crig(T) + crtAT) does not 
influence the results presented in the next section. 


To calculate the thermodynamic properties of the gluon sys¬ 
tem within the NAQPM we follow the procedure described in 
1 2511 to avoid thermodynamic inconsistencies. The energy den¬ 
sity is the sum of two contributions, coming from the trans¬ 
verse and longitudinal degrees of freedom respectively, given 
hy iji ^ \IT, kB^l): 


e — etr + eio 


where 




with 



exp[-——] 


gf,i aj{k, m) 
^ - 1 


i — tr, lo 


(8) 


(9) 


gfjr = 2-8 ; gfjo(k) = 8 ■ (1 - ^7^) (10) 

and 


Ni = NAni, a-) 


1 

cr 
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Counting the microstates we have included the prefactor 
and Ntr, Nh are normalization factors for the truncated 

normal distributions describing the gluon mass fluctuations for 
the transverse and longitudinal degrees of freedom respectively. 



In the counting of the number of degrees of freedom, we intro- 

k^ 

duce a momentum depending factor gfhik) = 8(1- : -r) 

k-^ +m‘^ 

for the longitudinal component which takes into account the 
fact that the longitudinal degrees of freedom vanish for large 
momenta. For the transverse degrees of freedom we use gf^tr - 
2 ■ 8 = 16. In both cases the 8 corresponds to the eight colored 
gluons. 

Setting k - T ■ X transverse and longitudinal energy con¬ 
tributions respectively become: 
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Substituting x - j ■ sinh(t) and using known properties of 
the Bessel functions, the integration over x is straightforward, 
leading to the expressions: 
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.{jlf.Ktijl) ( 11 ) 


with Ky the modihed Bessel functions 112811 . The mass integral is 
performed numerically. Having calculated the energy density, 
the pressure is obtained ll25n integrating the thermodynamically 
dP 

consistent relation e - T - P: 

dT 


P 

T 



( 12 ) 


The integral in Eq. ([12]) is also performed numerically. The final 
result for the energy density as well as the pressure, depends of 
course on the parameters fJ.tr{T), Pio{T) and cr(7’). 


4. Numerical results 

Using equations (18111112b we calculate the energy density 
and pressure for the gluon system at hnite temperature. For 
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temperatures T jTc ^ 1.85 we use the lattice data presented in 
||2^ for the transverse and longitudinal mean masses employing 
linear interpolation to determine their values in between. Be¬ 
yond TI Tc ~ 1.85 longitudinal and transverse masses approach 
each other and their common behaviour in this high temperature 
regime becomes linear following the relation; 

= 1.18(7’/?’,)+ 0.48 (13) 

obtained by a linear fit to the lattice data ll2^ in the region 
r > 1.5, with Tc ~ 270 MeV. The variance cr is determined by 
a Monte-Carlo search such that the lattice data for the energy 
density and pressure are fitted within 10 accuracy 0. Since 
there are two sets of lattice results for the gluon energy density 
density and pressure given in ||3l and |@] we have fitted both of 
them. In ||3l the equations of state are calculated for a relatively 
small temperature range above Tc while in 0 the results extend 
to the region of very high temperatures. Close to the critical 
region there is a small deviation between the two lattice calcu¬ 
lations probably due to finite size effects. Within our approach 
this deviation is reflected in the values of cr{T) necessary to fit 
the corresponding lattice results. 

In order to guide the search for the cr{T) profile so that an 
approximate continuous function emerges, we use the follow¬ 
ing strategy: Assume that we want to determine the parameter 
cr of our model for an (ordered) sequence of temperatures T,, 
i - \,2,..,N with Ti+i > T,. Furthermore, assume that we 
have located the optimal value cr(7’,) for the temperature T,. To 
estimate the value cr(7’,+i) at the subsequent temperature Ti+i 
we explore a region centered at cr(7’,) and extending up to 50% 
around it in each direction. This procedure turns out to con¬ 
verge surprisingly fast to the optimal value. This is crucial since 
the numerical calculations of the integrals in Eqs. (fTTT i are per¬ 
formed with high accuracy consuming CPU time. The results 
of our numerical analysis concerning the equation of state of 
the gluonic system are presented in Fig. 1. 

We observe an excellent ^reement between the lattice re¬ 
sults given in references [Sl lOfl and the calculations using the 
NAQPM introduced in the present work. Note that the agree¬ 
ment is maintained even in the high-temperature regime where 
the values approach from below the Stefan-Boltzmann limit, 
indicated by the dotted line in Fig. l(a,b). In Fig. 2(a,b) we 
show the dependence of the parameters nio (a) and cr (b) 
on the temperature T as obtained by the Monte-Carlo optimiza¬ 
tion procedure described above. In Fig. 2a we also display with 
points the lattice data. We emphasize that the linear dependence 
of the mass on T for T jTc > 1.9, as demonstrated in Fig. 2a, is 
implied by the lattice calculations in ll2^ . Notice that a similar 
functional form is obtained in high -T perturbative QCD 1^ . 
The red curve in Fig. 2b is obtained fitting the results of refer¬ 
ence ||3l while the black curve is the NAQPM result for the the 
description of the corresponding lattice data in reference 
The small shift in cr{T) is needed to capture the deviation be¬ 
tween the two lattice results. Furthermore, Figs. 2(a) and 2(b) 




Figure 1: (Color online) (a) The energy density scaled by T'^ of the gluon sys¬ 
tem as obtained by the NAQPM (red and black lines) and the lattice results 
for the same quantity (red asterisks and black circles), (b) The pressure (also 
scaled by T‘^) calculated with the NAQPM (shown with red and black lines) 
and the corresponding lattice results. The plots display the region T > Tc and 
t denotes the reduced temperature (T — Tc)/Tc. The lattice results are from H 
(red crosses) and d (black triangles). In both plots the dotted line indicates the 
Stefan-Boltzmann limit. 


^In fact the results are quite robust. Changing the convergence criterion 
from lO”"^ to 10“^ leads to a change in the second decimal digit for the param¬ 
eter (T. 
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indicate that in the high temperature regime the variance cr is 
proportional to the corresponding common mean mass, as both 
parameters become linear in T. This is a sign for the presence 
of a strong correlation between // and cr and indicates that the 
mass distribution becomes in fact mono-parametric. 




Figure 2: (Color online) (a) The ratios ^ and ^ as a function of the tem¬ 
perature T. With points we show the lattice data taken from |2^ . (b) The 
temperature dependence of the variance cr scaled by the critical temperature 
Tc- The plots display the region T >Tc and t denotes the reduced temperature 
{T - Tc)ITc- The black and red lines are obtained through a fit to the results of 
Ref. (i and Ref. fl respectively. 

A final comment is in order here. As seen in Fig. 2b, when 
approaching the critical point an increase of the variance cr is 
observed which does not appear in the plot of the masses shown 
in Fig. 2a. Thus, within NAQPM criticality is reflected in the 
temperature dependence of cr which defines naturally a measure 
for the mass fluctuations. This allows for an interpretation of 
the gluon deconfining transition as a change in the gluon mass 
spectrum. 

5. Concluding remarks 

We have introduced a quasi-particle model for the thermo¬ 
dynamical description of the gluon plasma which takes non- 
abelian characteristics into account. In particular, the mi¬ 
crostates building up the statistical ensemble for the gluon are 


determined from non-linear plane wave solutions of the asso¬ 
ciated classical dynamics, containing a free mass parameter. 
They correspond to quasi-particles with a continuously vary¬ 
ing mass having a quasi-Gaussian distribution with tempera¬ 
ture dependent mean mass and variance. To bridge the gap with 
the Lattice we use different mean masses for transverse 
and longitudinal gluonic degrees of freedom keeping a com¬ 
mon variance. Tuning appropriately the variance cr, which is 
the only free parameter in our model, we reproduce with high 
accuracy the results of the most recent Lattice calculations for 
the gluon plasma equation of state. At the same time we avoid 
singular behaviour of the mean mass close to the critical point, 
a feature which is common in QPMs but is in contradiction to 
Lattice results. In our approach the traces of the transition to the 
gluon plasma phase are imprinted on the rapid increase of the 
gluon mass variance in the neighbourhood of the critical point. 
It would be interesting to look for a similar behaviour in the 
phase diagram when matter degrees of freedom are included. 
However, this is left for future investigations. 
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